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We examine the time evolution of a two level ion interacting with a light field in harmonic 
oscillator trap and in a trap with anharmonicities. The anharmonicities of the trap are quantified 
in terms of the deformation parameter r characterizing the g-analog of the harmonic oscillator trap. 
Initially the ion is prepared in a Schrodinger cat state. The entanglement of the center of mass 
motional states and the internal degrees of freedom of the ion results in characteristic collapse and 
revival pattern. We calculate numerically the population inversion I{t), quasi-probabilities Q{t), 
and partial mutual quantum entropy S{P), for the system as a function of time. Interestingly, small 
deformations of the trap enhance the contrast between population inversion collapse and revival 
peaks as compared to the zero deformation case. For l3 — 3 and 4,( /3 determines the average number 
, of trap quanta linked to center of mass motion) the best collapse and revival sequence is obtained for 

■ T = 0.0047 and r = 0.004 respectively. For large values of r decoherence sets in accompanied by loss 

of amplitude of population inversion and for r ~ 0.1 the collapse and revival phenomenon disappear. 
Each collapse or revival of population inversion is characterized by a peak in S{P) versus t plot. 
During the transition from collapse to revival and vice- versa we have minimum mutual entropy value 
CN ; that is S{P) = 0. Successive revival peaks show a lowering of the local maximum point indicating a 

^ . dissipative irreversible change in the ionic state. Improved definition of collapse and revival pattern 

as the anharminicity of the trapping potential increases is also refiected in the Quasi- probability 
versus t plots. 



o 
o 

^ ' I. INTRODUCTION 



Two-level ions trapped in a harmonic or anharmonic trap have turned into an important tool for understanding the 
' time evolution of non classical states. Experimentally harmonic oscillator traps have been realized and various types 
^ ' of non-classical states of ions in trapped systems constructed [|l|,pl with sufficient control over amplitudes, relative 
O^. phases etc. of the component states. It is now possible to detect | J experimentally constructed non-classical states 
' and measure the statistics of the quantum motion of the center of mass of the ion. The non-classical states constructed 
and detected include the Fock states, thermal states, Shrodinger cat states and squeezed states. The possibility of 
manipulating the entangled states of particle systems provides the basis for applications like quantum computations 
^ , and quantum communications. On the other hand, these experimental advances have paved a way to an understanding 
- - ■ of some of the basic aspects of quantum theory. A controlled manipulation of these states requires precise control of 
the Hamiltonian of the system. The form of the trap potential is usually controlled by the geometry of the system, 
the laser field intensities and relative phases. As such anharmonic traps with different potential shapes should also be 
viable experimentally. The measurement of dynamical behavior of non-classical states in harmonic and anharmonic 
traps and a study of anharmonicity related new features makes an interesting study. Harmonically trapped and 
laser driven ion system has been studied in many contexts such as nonlinear coherent states vibrionic Jaynes- 

Cummings interactions j^, nonlinear Jaynes-Cummings interactions and generation of even and odd coherent 

states Q . The trapped ion-laser system is in some sense an experimental realization of Jaynes Cummings model , 
just like the micromaser cavity experiments. A study of the dynamics of two-level atom coupled to the q-analog of a 
single mode of the bosonic cavity field indicates that q deformation of Heisenberg algebra may correspond to some 
effective nonlinear interaction of the cavity mode. In our earlier work [p^ on ions in harmonic and anharmonic traps, 
we examined the collapses and revivals of the population inversion and quasi-probabilities as a function of time for the 
initial state of the system having ion in it's ground state and the center of mass motion described by a coherent state. 
This study brought out some very interesting new features in the dynamics of the system. In the present work, we 
consider a two- level ion initially in a Schrodinger cat state and examine the collapse and revival of population inversion 
and quasi-probabilities with a focus on the coherence of the system and anharmonicities of the trap. We also examine 
the time evolution of partial mutual entropy of the system to understand the improvement in collapse and revival 
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pattern of the system for characteristic values of parameter quantifying the trap anharmonicities. The motivation for 
studying these states is provided by the possible use in quantum computation and quantum communication where 
entanglement of ionic internal states with center of mass motion may be used to encode information ]r^-[l^ . In Ref. 
|p2| as well as in this work, the anharmonicity of the trap potential is modeled through a g-analog harmonic oscillator 
trap. Since Macfarlane and Bicdcnharn [|l8| discussed g-analog of harmonic oscillator, a lot of work using these 
in various fields of physics has been done. In our efforts to understand the physical nature of g-deformations in the 
context of pairing effect in nuclei [p^-pl| , the g-deformation is found to simulate the nonlinearities of the problem 
caused by unaccounted part of the residual interaction. In the present context, with the parameter q being used to 
quantify the anharmonicity of the trap potential, the underlying energy spectrum of vibrational states of ion as well 
as the laser ion interaction contain dependence on the number operator for vibrational quanta and the parameter q. 
For small values of g, the g-analog harmonic oscillator is essentially an anharmonic oscillator with anharmonicities 
. Nonlinear couplings of vibrational modes of ion have been studied theoretically p3|~p5| with the nonlinear effects 
modeled in the laser-ion interaction Hamiltonian through Hermitian operators with strong dependence on the number 
operator for vibrational quanta. The underlying trap in these studies is a harmonic oscillator one. We may point out 
that for ion in a harmonic oscillator trap and non-linear ion-laser coupling expressed through g— deformed operators 
with the g values used in our work, the system behavior continues to be very close to the case of ion in a harmonic 
trap and linear ion-laser couplings. 

We start with a description of the trapped ion-laser system Hamiltonian in section I for the cases of ion in a 
Harmonic Oscillator trap and ion in a g-analog harmonic oscillator trap. In section II the initial state and the 
equations governing it's time evolution are given. Population inversion and quasi-probability are defined in section 
HI. In section IV, we discuss the concept of partial mutual entropy for the system and the numerical results are 
analyzed in section V. 



II. THE HAMILTONIAN 



A. Harmonic Oscillator trap 



The system that consists of the two level ion moving in a harmonic oscillator trap potential and interacting with a 
classical single-mode light field of frequency uji is described by the Hamiltonian, 



H 



2 



F*a-) 



(1) 

uji , is the detuning 



where a and are the usual creation and destruction operators for the oscillator, A = t 

parameter and f2 is the Rabi frequency of the system. The pseudo-spin operators of the two-level ion, tr^ and az 
operate on state vectors \g) and |e) where indices g and e stand for the ground and excited states of the two level ion. 
The operator F stands for Q exp(iA:x) = exp[ie(a^ + a)], where k is the wave vector of the light field and x position 

quadrature of the center of mass. The parameter e — \J^^ is a function of the ratio of the recoil energy of the ion 

Ej. = and the characteristic trap quantum energy Ft = Tiaj — The second term in the Hamiltonian refers to 
the energy associated with internal degrees of freedom of the ion, whereas the third term is the interaction of the ion 
with the light field. By using a special form of Baker-Hausdorff Theorem the operator F may be written as a product 
of operators i.e.e we use the equality 



F — exp[ie(a^ + a)] = e'" 
For harmonic oscillator trap one has 

[a,a^] = l; N ~ N ^ ; N a - 
The physical processes implied by the various terms of the operator 



^a])g^eatg^ea 



aN 



F = 



exp 



— e 



E 



it" 



■E 



fc! 



(2) 



(3) 



(4) 



may be divided into three categories. The terms for n > k correspond to an increase in energy linked with the motional 
state of center of mass of the ion by {n — k) quanta. The terms with n < k represent destruction of {k — n) quanta 
of energy thus reducing the amount of energy linked with the center of mass motion. For (n = k), we have diagonal 
contributions. The contribution from a particular term containing operators a^^a'^is determined by the coefficient 



exp 



(^) 



n!fe! 
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B. g-analog harmonic oscillator trap 



Next we consider a single two level ion having ionic transition frequency uoa in a quantized g-analog quantum 
harmonic oscillator trap (g-deformed harmonic oscillator trap) interacting with a single mode travelling light field. 
The creation and annihilation operators for the trap quanta satisfy the following commutation relations, 



AA^ -qA^A^q-^ ■ NA^-A'^N^A^ ■ NA- AN = -A. 
The operators A and A^ act in a Hilbert space with basis vectors |n),rt = 0,l,2,..., given by, 



-|o) 



such that the number operator N \n) ^ n\n) . The vacuum state is A |0) =0. We define here [x\q as 

— q^^ 



(5) 



(6) 



(7) 



and the q- analog factorial [n\q\ is recursively defined by 

[0]ql = [l]g! — 1 and [n]ql — [n]q[n — l]ql. It is easily verified that 



A^\n) = + |n+l) 



A\n) = [n]^ \n-l) 



(8) 



and N is not equal to A^A. Analogous to the harmonic oscillator one may define the g-momentum and q-position 
coordinate 



^^r-^iA^~A) 



2muj 



{A^ + A). 



The g-analog harmonic oscillator Hamiltonian is given by 



Hqho = ^nw{AA^ +A^A) 



(9) 



(10) 



with eigenvalues 



En = -^(^([n- + l]g + Mq 



(11) 



We note that the trap states are not evenly spaced, the energy spacing being a function of deformation. Besides that 
as we move up in the number of vibrational quanta in the states the spacing between successive states increases. We 
choose q = , where r is a real valued parameter. For r — 0.0 the harmonic oscillator trap is recovered. As such the 
parameter r is a measure of the extent to which the harmonic oscillator potential trap is deformed. Using the Taylor 
expansion of the q - number \n]q in terms of powers of the energy eigenvalues may be rewritten for small r as 



E,-, = hu 



T 

24 



1 



6 



(12) 



Bonatsos et. Al 
looks like 



have shown that the potential giving a spectrum similar to that of Eq. (^) up to the order 



V{x) = 



120 



(13) 



which is an anharmonic oscillator with anharmonicities 

Using a nonlinear map given by Curtright and Zachos [E? | one can express Hqho in terms of the operators a and 
of the harmonic oscillator. To make this point clear we express the operators A and A'^ as 



A^afiN) 



At = /(AT) at 



(14) 
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where f{N) — ^-^^^^ ^ and N — a^a. We can also verify that 

/(iV)at = at/(iV + 1) ; f{N)a = af{N - 1) (15) 
The Hamntonian of Eq. can be rewritten as 

Hgho = hiv {f{N + If + f{Nf) (^atfl + (16) 

which can be interpreted as a harmonic oscillator Hamiltonian with a frequency ujq{N) that depends on the quantum 
number n of the state in question. 

The quadratures x and p are related to Pq and Xg through 



Pq=^^j'^ [[f{N) !{N + 1)) ^ - {f{N) + m + 1)) ^ 
X, = ^ ((/(iV) + fiN + D) - ifiN) - /(A. + D) ^) (17) 



C. Hamiltonian for an ion interacting with light in a g-analog harmonic oscillator trap 

The Hamiltonian for an ion interacting with light in a q-analog harmonic oscillator trap may now be written as 

Hq = i?iw(AAt + A^A) + ^hAa^ + ^nn{Fq(j+ + F^a-) 



where by analogy with Eq. |j we choose 
which in the limit q ^ 1 reduces to F, and can be expanded as 



= eV ' Je'^^'e''^ , (18) 



j-^) (»e)"At" ^ (ze)Mfc 



^.-V ^^l^^^^l^^^. (19) 



n=0 fc=0 



Various terms in the expansion of this operator represent processes which might result in transitions of the center of 
mass from a given motional state, in the g-analog trap, to another, while loosing or gaining energy. The coefficient of 

the operators A^^'A'' is again exp ''^'ri\k\ ' same as that in the corresponding term with operators a'''"a'^ 

in Eq. ^. Intuitively this is the correct way of representing the interaction of the ion and the laser in a g-analog 
harmonic oscillator trap. This form of the operator shows that the energy exchange of the center of mass motion of 
the ion occurs as the ion moves up or down in the trap. Since En+i — En is n dependent, it implies an n dependent 
entanglement of the center of mass motion and the internal degrees of freedom of the two level atom. Using Eq. |l^, 
we can rewrite Eq. ^ as 

^ - - (,)n+fe (e/(jV)at) " (6«/(iV))^ 

n=0 k=0 

Comparing with F we notice that the effective lamb Dicke parameter in a g-analog trap, for the loss and gain of 
motional state energy in an interaction process, is e/(iV) where N is the number operator. With the g-analog of 
Glauber coherent state defined as 



/expl"! n=o J[n]q\ 



oo „ 

a 



4 



the calculated value of q (a| f{N) \a) ^ is 1.0004, for the choice a = 4 , e = 0.05 and t = 0.003. 

It is important to recall at this point that the operator {A^ + A) is not a self adjoint operator. The vectors 

\m) are analytic vectors of the operators {A'^)'^A^ but not of the operators A'^(A'I')". One can easily verify that 

the matrix elements (m| e V ^ \^h^o '^n'=o — |w) Eire not well defined. On the other hand the matrix 

elements of the operator in Eq. |l^ are well defined as discussed in Ref. ||2^ . From the form of the operator in Eq. |l^ 
it is clear that the product _F)jCr "'"contains all the processes in which the atom is excited from the ground state to the 
excited state while at the same time some of the energy quanta in the center of mass motion are lost or gained. It is 
also important to note that Xq is not the true position coordinate but is related to x and p in some complex way. As 
such the operator exp[ie(A^ + A)] can not be used to represent the laser-ion interaction. 

The expression for the matrix element of the operator Fq between the states (m| and \n) for to < n is given by, 

N|! f-^^ kl{n - m + k)l[m ^ k]ql 



III. SCHRODINGER CAT STATE 



Initially the ion is prepared in a Schrodinger Cat State with equal probabilities of finding the ion in its ground state 
and in the excited state coupled to a coherent state describing the state of motion of the center of mass of the ion in 
the anharmonic trap. The experimental technique for producing trapped ion in a state with these characteristics has 
been described in many experimental papers and such states have been produced [^,^ . 

The initial state of the trapped ion in a g-analog harmonic oscillator trap can be expressed as 



\g,(3)+e"t' |e,-/3)„ 
^ft = 0)= " ^ ' (22) 



\g, (3)^ = E ^r= If' ' (23) 



where \g,P)q is given by 



expf "=0 ^N,! 

and (3 determines the average number of vibrational quanta linked to the center of mass motion. The relative phase 
between two possible internal states is determined by 0, in the present study we choose e*'^ = 1. 
The state of the system at a time i, 

■m m 

is a solution of the time dependent Schrodinger equation 



i?,*(t) = (25) 



In the state \g,m) (|e,m) ) the two level ion in it's ground (excited) state is coupled to number state |to) of the 
anharmonic oscillator. The probability amplitudes gmit) and emit) satisfy the following set of coupled equations 



.dgm{t) 
I 



\gm{t)H[ni + l]q + [m]q)- /\) + ]pY.en{t){g,m\'^ Fl\e,n) 

n 

'^^^ = lernit) {uj{[m + l]q + [m]q) + A) + if] ^ 5„(t) (e, m\ Fqa+ |g, n) (26) 

n 

In the limit g — > 1 we recover the system that gives the dynamics of two-level ion in a harmonic oscillator trap. By 
rescaling e, w, and t, the parameter can be eliminated from Eq. (|2q). 
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IV. POPULATION INVERSION AND QUASI-PROBABILITY 



The entanglement of motional degrees of freedom of the center of mass of the ion with it's internal degrees of 
freedom manifests itself in the well known collapse and revival of population inversion. The population inversion is 
defined as 

m=Pg{t)-P,{t), (27) 

that is the difference between the probability of finding the system in the ground state, Pg{t), and the probability of 
finding the system in the excited state, Pe{t). The time evolution of population inversion is examined for different 
values of trap anharmonicities for the cases when coherent state is characterized by parameter /? = 3, 4 . The time t 
is expressed in units of 

We have also calculated numerically and plotted as a function of ar and the Quasi-Probability function Q (a) 
defined as 

Q{a)^- {a\pc.m.{t)\a) , (28) 

TTq 

where pc.m.it) = trionp(t) is the density matrix reduced for degrees of freedom of center of mass motion. 



V. QUANTUM MUTUAL ENTROPY 



Dynamical behavior of the system can be better understood in terms of quantum mutual entropy. Furuichi and 
co-workers have applied the concept of quantum mutual entropy to dynamical change of state of the atom in 
Jaynes-Cummings model(JCM) The analogy between JCM and ion in a trap system allows us to extend the 

concept to understand the entanglement of the internal degrees of freedom of ion and it's motional degrees of freedom. 
We next outline the calculation of quantum mutual entropy for the system at hand prepared initially in state given 
by Eq. ( p2[ ) . It is found that the quantum mutual entropy can be decomposed into a part determined by the relative 
populations of the system and another that depends on the coherences. The partial quantum mutual entropy used 
by us in our numerical study is essentially a measure of mutual entropy due to populations developed in the system 
with the passage of time. 

With the initial state of the trapped ion expressed as in (p^). The density operator for the internal states of the 
ion at i = is 

P*on = Ai|g)(.g|+A2|e)(e|. (29) 

The state of motion of the ionic center of mass at t = is represented by a coherent state for ion in ground state 
and by coherent state |— /?)^ for ion in excited state, where 

E^l-), ■ (30) 



exp„l|/?ri "=o 



The state of the system at a time t is given by Eq. (24) 



The time evolution of the internal state of the ion can also be represented by a continuous mapping 

A^P^on = tr,,ra. {UtPoUt) (31) 

where po is the initial state of the system, Ut — exp ( ) a unitary operator and the mapping A"" maps the 
initial internal state pion to the current internal state of the ion at time t. We represent the time evolution of initial 
constituent ionic states \g) {g\ and |e) (e| through the maps 

AT \g) (5l = tr^.m. {Ut |.g,/3) {g,(3\ Ut) = P^'\t) \g) {g\ + Pi'\t) \e) (e| (32) 
+C«(i)|5)(e|+cW*(t)|e)(g| (33) 



and 



|e) (e| = tr,.„, {Ut |e, -(3) (e, -(3\ Ut) . (34) 
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The correlation between the initial ionic state and its current state containing information about the time evolution 
of each constituent ionic states is given in terms of the quantum mutual entropy defined as 



t Hion J 1 



(35) 



1=1 



where the relative quantum entropy 



S (Ar N) (*l , Ar p^on) = tr (Ar K> h (Ar \^) 

-log(Arfton)] 



(36) 



and i = 1 (2) for ion in ground (excited) state. It measures the relative overlap of the constituent ionic state \i) with 
the ionic state of the system at a time t for the cases a) where the initial ionic state of the system is state \i) and b) 
the initial ionic state of the system is given by Eq. (29). The quantum relative entropy is a more general concept than 
the Von Neumann Entropy being a measure that extends to mixed states. As defined in Eq. {Bw, it is a measure of 
how far each entangled state is from its component disentangled state as a result of time evolution. The quantum 
mutual entropy as defined by Eq. ^ is obtained by summing up the time dependent quantum relative entropies with 
respect to ionic states that constitute the initial cat state and measures the propagation of information content of the 
initial state. 

Working in the space spanned by basis vectors {|itn)} = l^j^^-) 7 ,r7i = 0, 1, 00, we find from Eq. ( p^ the 

probabilities of finding the ion in the ground state and the excited state (the populations) to be respectively 



(37) 



and 



P,{t) = (e \tr,onPit)\ e)=J2 K{t)\' 



(38) 



Besides that the system at a time t is characterized by the following non-diagonal matrix element of the operator 
trionPit) (coherence), 



Cge = {9\trionp{t)\e) =^5*(t)e„(t) 



(39) 



We also define the populations and coherences for the system prepared in initial state \g, (3) ^ {i = 1) and the initial 
state with the ion in the excited state |e, —f3)^ {i — 2). The state of the system at a time t is now given by 



(40) 



with the populations defined as 



pW(i) = 



tnonP^'Ht) 



(41) 



and 



pW(t) = (e|tr„„pW(i)|e)=Ehn^W 

n 

The coherences for these initial states are given by 

Ci'^ it) - (.9 |t^.o„P« it) I e) = E {9^'^ (^)) * it) 



(42) 



(43) 



We can verify that in terms of the populations and the coherences defined above the relative quantum entropy, 



A'f^Pion) (Eq. ^) for the system starting in initial state given by Eq. (|2|) is 
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pW^logf^j+pWlog 



Peit) 



-2 * Re 



(44) 



Substituting Eq. (M3) in the definition of quantum mutual entropy for the system, we have 



1=1.2 



(0/ 



-2* Re 
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(45) 



We can spUt Sm in to two distinct parts, one depending on populations and the other on off diagonal coherences that 



S{P) + S{C) 



(46) 



where 



S{P) = ^ A(i) 



1=1,2 



^'i"«)ioJ^)+pi"(mog(||> 



(47) 



We have calculated numerically S{P) the partial mutual quantum entropy for the cases /? = 3,4 and used the S{P) 
peaks to pinpoint the t values for which the system is the most correlated and the t values where it shows the least 
correlation. Essentially S{P) is used as a measure of entanglement of the system. 



VI. RESULTS 

The population inversion as a function of rescaled time parameter t{^) for two level ion interacting with light 
field in a harmonic oscillator trap and a g-deormed oscillator trap for the cases /3 = 3, 4 is displayed in Figures. 1 
and 3 for three different values of deformation parameter t {q = exp(T) with r real). The set of parameters used in 
the numerical calculation is cJ = ^ = 50, e = n ^ 0.05, and A = = —50. The maximum value of n in Eq. ( [2^ ) 
is restricted to rt = 32. For both /? values the collapse and revival pattern was found to improve in definition with 
increase in the anharmonicity of the trap potential, measured by parameter t. In Figures, la and 2a, we display 
the collapse and revival for ion in a harmonic oscillator trap that is r = 0.0. Both for f3 — 3 and 4 first revival and 
second revival are easily identified. For (3 = 3 the second collapse, which is barely identifiable is succeeded by a region 
where the inversion pattern becomes diffuse, whereas for /3 = 4 the second collapse is followed by another revival 
which is quite wide. For (3 = 3 and 4, the best collapse and revival sequence is obtained for r = 0.0047 and r = 0.004 
respectively as shown in Figures, lb and 2b. We can identify four revivals in each case. The collapses and revivals 
for P = A are remarkably sharp. Further increase in r results in a loss of contrast between successive revival peaks. 
Figures. Ic and 2c indicate a rapid loss of phase relationship essential for generating observable collapses and revivals 
of population inversion between components of the total wave function as the anharmonicity of the trap potential is 
increased. We have also found that for a large value of r, different for each /3 value, not only the collapses and revivals 
of population inversion vanish but also the time dependence of population inversion disappears in all cases. 

In figs. (3-4) we plot the partial quantum mutual entropy S{P) as a function of time for initial states with /3 = 3, 
and 4. A careful examination shows that each collapse as well as revival of population inversion is characterized by a 
peak in S{P) versus t plot. During the transition from collapse to revival and vice- versa we have minimum mutual 
entropy value that is S{P) = 0. In the q-deormed oscillator trap for the cases /3 = 3, r = 0.0047 and /3 = 4, r = 0.004 
we notice relatively higher S{P) peaks and well defined transition periods(collapse <;4> revival) in comparison with 
the harmonic oscillator trap(T = 0.0). The regions without continuous minimum S{P) do not show well defined 
collapses or revivals in the population inversion plots. In Table I, the t values for successive peaks in S{P) for the 
cases (3 = 3,T = 0.0047 and (3 = 4:,t = 0.004 are displayed. It is evident from figures 3 and 4 that S{P) decreases with 
the Rabi oscillations. Table I demonstrates that the successive revival peaks show a lowering of the local maximum 
point indicating a dissipative irreversible change in the ionic state. 
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Table I. The values of t for successive peaks in S{P) for the cases 
/3 = 3, r = 0.0047, /? = 4, t = 0.0 and /? = 4, t = 0.004. 



^ = 3,T = 0.0047 


/3 = 4,r = 0.0 


/? = 4,r = 0.004 


t 


S(P) 


Peak type 


t 


S{P) 


Peak type 


t 


SiP) 


Peak type 


0.0 


0.693 


Initial state 


0.0 


0.693 


Initial state 


0.0 


0.693 


Initial state 


58.6 


0.628 


revival 


85.8 


0.333 


revival 


67.8 


0.731 


revival 


114.0 


0.314 


collapse 


171.4 


0.143 


collapse 


133.2 


0.497 


collapse 


175.6 


0.285 


revival 


266.8 


0.106 


revival 


201.2 


0.512 


revival 


234.0 


0.220 


collapse 


388.2 


0.085 




266.6 


0.359 


collapse 


295.6 


0.181 


revival 


447.6 


0.076 




336.8 


0.397 


revival 


360.0 


0.234 


collapse 








404.8 


0.334 


collapse 


415.4 


0.175 


revival 








472.6 


0.324 


revival 



Next we plot in Fig. 5 the Quasi-probabilities for /? = 4 at t = 0.0,85.8,171.4,266.8,388.2 and 447.6 to visualize 
the time evolution of ionic center of mass in a harmonic oscillator trap(r = 0.0). As expected, the coherent states 
\g,l3) and |e,— /?) , each splits into two and move in a sense opposite to each other in phase space. At t — 85.8 
there are two somewhat distorted coherent states characterized by almost equal average position coordinate(ar — 0), 
and momenta with opposite signs. This point corresponds to a revival peak in the population inversion plot and the 
first peak S(P) versus t plot( Fig. 4a). A second revival similarly occurs aX t = 266.8. The collapse at t = 171.4 
corresponds to a situation similar to the initial state. At 388.2 and 447.6 the quasi-probability is spread out in the 
phase space. For r = 0.004 quasi-probability contour plots a,tt = 0.0, 67.8, 133.2, 201.2, 266.6, and 336.8 are presented 
in Fig. 6. Characteristic spreading of the Quasi-probabilities in the phase space is the result of non-linearities of the 
trap, however the initial state configuration is still recovered signalling collapse and revival phenomenon. 

VII. CONCLUSION 

Entangled quantum systems present correlations having no classical counterparts and may be exploited to store 
and transfer information in ways otherwise difficult or impossible. In this work we have examined the effects of trap 
non-linearity on Schrodinger cat state of a single trapped ion. This particular quantum system may serve in future as 
a basic element in quantum computation and quantum communication p"3|~p^ or a sensitive detector of decoherence 
due to interaction with environment. The entanglement of ionic internal states with center of mass motion may be 
used to encode information in such applications. We have found an enhancement in the contrast between the collapse 
and revival peaks as the anharmonicity of the trapping potential increases. With further increase in non-linearity 
not only the collapses and revivals disappear but also the time dependence of population inversion vanishes. This 
extremely interesting feature is manifested in the Population inversion plots, the Quasi- probability plots as well as 
the partial mutual entropy plots. The effect is markedly visible in f3 — 4 case for the set of parameters used to 
characterize the ion-laser in trap system. The best collapse and revival sequences are obtained for r = 0.0047 and 
r = 0.004, for (3 = 3 and 4 respectively. The reason for this effect is to be found in the deformation dependence 
of trap energies. The characteristic Rabi Frequency, in addition to the ion-laser interaction now depends on energy 
separation between trap levels. It can be easily verified that in the lowest approximation, with uj ^ fl and driving 

field tuned to A = — cjjthe Rabi frequency is given by /i(n) = (^(cosh(2(n + 1)t) — 1))^ + il'^e'^[n + l]q which in 

the limit r ^ (for q = e^) reduces to ^{n) = ^J^'^e'^[n + 1] as expected. For small values of parameter t, the factor 

(^(cosh(2(n -|- l)r) — 1))^, coming from non-equidistant trap energy spectrum, varies approximately as (n -|- 1)** and 
dominates the scene resulting in a high contrast collapse and revival peaks in population inversion. For large values of 
T decoherence sets in accompanied by loss of amplitude of population inversion and for r 0.1 the collapse and revival 
phenomenon disappear. We use partial quantum mutual entropy as an indicator of time evolution of entanglement. It 
is obtained by summing up the time dependent quantum relative entropies with respect to ionic states that constitute 
the initial cat state and retaining the part that depends on populations. Quantum relative entropy is a measure of 
how far each entangled state moves from its component disentangled state as a result of time evolution. Looking 
together at the population inversion and partial mutual entropy S{P) plots we verify that the onset of collapse and 
revival is characterized by S{P) = 0. A well defined collapse and revival pattern is characterized by wide S{P) = 
regions. Every collapse or revival corresponds to a peaked S{P). Successive revival peaks show a lowering of the local 
maximum point which is an indicator of a dissipative irreversible change in the ionic state same being true for S{P) 
peaks that indicate collapse of the population inversion. 
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g. 1. Population Inversion I{t) versus f for /3 = 3 and r = 0.0, 0.0047, 0.008. 
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Fig. 2. Partial mutual quantum entropy S{P) versus t for/3 = 3 and r = 0.0, 0.0047, 0.008. 



12 




13 



0.4 
0.0 
0.4 
0.0 
0.4 
0.0 



1 ' 

1 p=4, X=0.0 (a) 

j ill 






x=0.00 

1 1 1 

1 1 1 1 


4 


(b) 

JiilwiwJijilliL 




\ 

x=0.00 


18 









200 



400 



Fig. 4. Partial mutual quantum entropy S{P) versus t for /3 = 4 and r = 0.0, 0.004, 0.008. 
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Fig. 5. Quasi-probability plots for /3 = 4 and r = 0.0 att = 0.0,85.8,171.4,266.8,388.2 and 447.6. 
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Fig. 6. Quasi- probability plots for /3 = 4 and t = 0.004 att = 0.0, 67.8, 133.2, 201.2, 266.6, and 336.8. 
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